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Abstract. As it has been proven, the determination of general one-dimensional Schrodinger 
Hamiltonians having third-order differential ladder operators requires to solve the Painleve 
IV equation. In this work, it will be shown that some specific subsets of the higher-order 
supersymmetric partners of the harmonic oscillator possess third-order differential ladder 
operators. This allows us to introduce a simple technique for generating solutions of the 
Painleve IV equation. Finally, we classify these solutions into three relevant hierarchies. 
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1 Introduction 

Nowadays there is a growing interest in the study of nonlinear phenomena and their corre- 
sponding description. This motivates to look for the different relations which can be established 
between a given subject and nonlinear differential equations pp. Specifically, for supersymmetric 
quantum mechanics (SUSY QM), the standard connection involves the Riccati equation [21 [3], 
which is the simplest nonlinear first-order differential equation, naturally associated to general 
Schrodinger eigenproblems. Moreover, for particular potentials there are other links, e.g., the 
SUSY partners of the free particle are connected with solutions of the KdV equation [4j. Is there 
something similar for potentials different from the free particle? 

In this paper we are going to explore further the established relation between the SUSY 
partners of the harmonic oscillator and some analytic solutions of the Painleve IV equation {Pjy)- 
This has been studied widely both in the context of dressing chains [51 El E] and in the 
framework of SUSY QM [9l [TO] [TTj 112] . The key point is the following: the determination 
of general Schrodinger Hamiltonians having third-order differential ladder operators requires 
to find solutions of Pjy. At algebraic level, this means that the corresponding systems are 
characterized by second-order polynomial deformations of the Heisenberg-Weyl algebra, also 
called polynomial Heisenberg algebras (PHA) [IT]. It is interesting to note that some generalized 
quantum mechanical systems have been as well suggested on the basis of polynomial quantum 
algebras, leading thus to a g-analogue of Pjy [8]. 

On the other hand, if one wishes to obtain solutions of Pjy, the mechanism works in the 
opposite way: first one looks for a system ruled by third-order differential ladder operators; 
then the corresponding solutions of Pjy can be identified. It is worth to note that the first-order 
SUSY partners of the harmonic oscillator have associated natural differential ladder operators of 

*This paper is a contribution to the Proceedings of the Workshop "Supersymmetric Quantum Me- 
chanics and Spectral Design" (July 18-30, 2010, Benasque, Spain). The full collection is available at 
|http://www.emis.de/journals/SIGMA/SUSYQM2010.html1 
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third order, so that a family of solutions of F\y can be easily obtained through this approach. Up 
to our knowledge, Flaschka (1980) P3] was the first people who realized the connection between 
second-order PHA (called commutator representation in this work) and P\y equation (see also 
the article of Ablowitz, Ramani and Segur (1980) [H]). Later, Veselov and Shabat (1993) [5], 
Dubov, Eleonsky and Kulagin (1994) [6J and Adler (1994) [7J connected both subjects with first- 
order SUSY QM. This relation has been further explored in the higher-order case by Andrianov, 
Cannata, Ioffe and Nishnianidze (2000) [9], Fernandez, Negro and Nieto (2004) [ID], Carballo, 
Fernandez, Negro and Nieto (2004) p], Mateo and Negro (2008) [12], Clarkson et al. [El EE], 
Marquette [TT] [18], among others. 

The outline of this work is the following. In Sections [2HH] we present the required background 
theory, studied earlier by different authors, while Sections [6j [7] contain the main results of this 
paper. Thus, in the next Section we present a short overview of SUSY QM. The polynomial 
deformations of the Heisenberg-Weyl algebra will be studied at Section [3l In Section 0] we 
will address the second-order PHA, the determination of the general systems having third-order 
differential ladder operators and their connection with Painleve IV equation, while in Section [5] 
the SUSY partner potentials of the harmonic oscillator will be analyzed. In Section [6J we will 
formulate a theorem with the requirements that these SUSY partners have to obey in order to 
generate solutions of Prv, and in Section [7] we shall explore three relevant solution hierarchies. 
Our conclusions will be presented at Section [HJ 



2 Supersymmetric quantum mechanics 

Let us consider the following standard intertwining relations [2[ [3j HI EDJ (2TJ [22] 



HiAf = Af Hi-!, Hi-iAj = At Hi 



A ± - — 



d , . 
T-r- + oii{x,ei) 
ax 



k. 



where 



Hi 



1 d 2 



+ Vi(x), 



i = 0, . . . , k. 



2dx 2 

Hence, the following equations have to be satisfied 

c4(x, €i) + of (x, €i) = 2[Vi-i(x) - €i], 
Vi(x) = 14-1 (re) - a-(x,ej). 



(2.1) 
(2.2) 



The previous expressions imply that, once it is found a solution ai(x,€i) of the Riccati 
equation (|2.ip associated to the potential Vi-i(x) and the factorization energy the new 
potential Vi(x) is completely determined by equation (|2.2[) . The key point of this procedure is 
to realize that the Riccati solution oti(x, e^) of the i-th equation can be algebraically determined 
of two Riccati solutions aj_i(x, ei_i), ai~i(x, ei) of the (i — l)-th equation in the way: 

/ \ / \ 2(ej_i - €j) 

aj_i(x,ej_i) - ai-i{x,ei) 

By iterating down this equation, it turns out that ai(x, ei) is determined either by i solutions of 
the initial Riccati equation 



a[(x,ej) + a\{x,ej) = 2[V (x) - ej], j = 1, ... ,i, 
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or by i solutions Uj of the corresponding Schrodinger equation 

H Uj = --u" + V (x)uj = ejUj, j = 1, . . . , i, (2.3) 

where a±(x,€j) = u'j/uj. 

Now, the k-th order supersymmetric quantum mechanics realizes the standard SUSY algebra 
with two generators 

[Qi,H ss ]=0, {Q l ,Q J } = ^H ss , i,j = 1,2, 

in the way 



Q + 


-(? 


H d 


-ei)- 



H ss = {Q^,Q + }= ( B *f* =(H d -ei)-(H d -£ fc ), (2-4) 

where 

i_i \ 

H d -e,= ^ Q ffo _ e J, , = 1,...,*. 

Notice that, in this formalism, Hq and H k are the intertwined initial and final Schrodinger 
Hamiltonians respectively, which means that 

H k B+ = B+H , H B^ = B^Hk, (2.5) 

with 

B t = A~l ■ ■ ■ Af, B k = A 1 ■■■A k , 

being k-th. order differential intertwining operators. The initial and final potentials Vq, V k are 
related by 



V k {x) = V (x) - ^) = Vo(x) - {\n[W( Ul , . . . , u k )}}", 



i=l 



where W(u\, . . . ,u k ) is the Wronskian of the k Schrodinger seed solutions Uj, which satisfy 
equation (j2.3|) and are chosen to implement the A:-th order SUSY transformation. 

3 Polynomial Heisenberg algebras: algebraic properties 

The PHA of m-th ordei0 are deformations of the Heisenberg- Weyl algebra [21 El IS] : 
[H, L ± ] = ±L ± , [L-,L + ] = Q rn+1 (H + 1) - Q m+ i(H) = P m (H), 

m+l 

Q m +i( H ) = L+L ~ = HiH-£i), 

i=l 

where the Hamiltonian describing those systems has the standard Schrodinger form 



1 In this work we will use the terminology of [TT], although in [23] these systems are called of (m + l)-th order. 
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■ £j + I - 1 

u J; + l ns, + i "•• •£. + ! a£ ; + 1 ■■• aSj + 1 ■■• ms.+i 
mSl me, a£l m £j 

(a) (b) 

Figure 1. Possible spectra for systems characterized by m-th order PHA, s < m + 1. In case (a), the 
spectrum consists of s infinite physical ladders, while in case (b) the j-th ladder has finite length since 
equation (|3.ip is satisfied. 



Q m +i{x) is a (m + l)-th order polynomial in x which implies that P m (x) is a polynomial of 
order m in x, and in this work we assume that L~ = (L + )^ . 

The corresponding spectrum of H, Sp(-ff), is determined by the s physical eigenstates of H 
belonging as well to the kernel of L~ : 

L~ips t = 0, Hip £i = £i^£ tl i = 1, . . . , s, 

where s is an integer taking a value between 1 and m + 1. s independent ladders can be 
constructed departing from these physical extremal states by acting repeatedly the operator L + 
on each ip£.. In general, these ladders are of infinite length. An illustration of this situation is 
shown in Fig. [T]a. 

On the other hand, for a certain extremal state tpg. it could happen that 

(L+y-^g-jkO, (L+) l i> £j =0, (3.1) 

for some integer L In this case it turns out that one of the roots different from the physical ones 
{Si, . . . , £ s } acquires the form 

£i = £j + l, i £ {s + 1, ...,m + 1}. 

Thus, the j-th ladder, which departs from £j, will end at £j + 1 — 1, i.e., it has finite length. An 
illustration of this situation is shown in Fig. Q]&. 



4 Second-order polynomial Heisenberg algebras 

By making m = 2 we get now the second-order PHA for which O [5l El [24"1 [25] : 

Q 3 (H) = {H- Sx) (H - £ 2 ) (H - £ 3 ) , 

P 2 (H) = 3H 2 + [3 - 2(£i + £ 2 + £ 3 )]H + l-{£ 1 + £ 2 + £ 3 ) + £ x £ 2 + £ x £ 3 + £ 2 £ 3 . 



(4.1) 



From the purely algebraic point of view, the systems ruled by second-order PHA could have up 
to three independent physical ladders, each one starting from £\, £ 2 and £ 3 . 

Now, let us take a look at the differential representation of the second-order PHA. Suppose 
that L + is a third-order differential ladder operator, chosen by simplicity in the following way: 



T + T + 



1 

71 



d 2 , . d . . ' 
— +9 (x)- + h(x) 



These operators satisfy the following intertwining relationships: 



HLt = L+(H a + l), 



L+H 



[H,L- 
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where H a is an intermediate auxiliary Schrodinger Hamiltonian. By using the standard equations 
for the first and second-order SUSY QM, it turns out that 

-f' + f = 2(V -St), 
V a = V + f'-l = V + g', 

Jl / „/ \ 2 



9" I 9' Y i , 9 2 , (£2 ~ S 3 ) 2 



25 \29 



r 



+ S 2 + S 3 -2 = 2V, 



q< q 2 

h = - y - + y --2V + S 2 + S z -2. 



By decoupling this system it is obtained: 



/ = x + 9 
h ■ 



2 , 9 
X + 2 



2 



2xg + a, 



rr x 9 9 1 

F = — + — +xq + Si , 

2 2 2 2 



(4.2) 
(4.3) 

(4.4) 



where 



2.9 



+ -g 3 + 4xg 2 + 2 (x 2 



Notice that this is the Painleve IV equation (F\v) with parameters a = S2 + S 3 — 2£i — 1, 
& = -2(£ 2 -£ 3 ) 2 . 

As can be seen, if one solution g(x) of Prv is obtained for certain values of Si, S2, S3, then the 
potential V(x) as well as the corresponding ladder operators are completely determined (see 
equations (|4.2p ~ (|4.4p ). Moreover, the three extremal states, some of which could have physical 
interpretation, are obtained from the following expressions 



ip £l tx exp ( — - 



, 9' 9 A 

,9' 9 , A 

ip£, oc 

1 2g 2 g 



gdx 



x l exp 



exp J 



9^ + 9_A 
29 2 g 



2.9 



9 ^ 
- + — 

2 5 



(4.5) 



dx 
dx 



where A = Si — S3. The corresponding physical ladders of our system are obtained departing 
from the extremal states with physical meaning. In this way we can determine the spectrum of 
the Hamiltonian H. 

On the other hand, if we have identified a system ruled by third-order differential ladder op- 
erators, it is possible to design a mechanism for obtaining solutions of the Painleve IV equation. 
The key point of this procedure is to identify the extremal states of our system; then, from the 
expression for the extremal state of equation (|4.5p it is straightforward to see that 



g(x) = -x - {ln[^ £l (x)}}' . 



Notice that, by making cyclic permutations of the indices of the initially assigned extremal states 
V"£i> ^£ 2 ) ^£31 we obtain three solutions of Prv with different parameters a, b. 
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5 k-th order SUSY partners of the harmonic oscillator 

To apply the k-th order SUSY transformations to the harmonic oscillator we need the general 
solution u(x,e) of the Schrodinger equation for Vq(x) = x 2 /2 and an arbitrary factorization 
energy e. A straightforward calculation [26] leads to: 



u[x, e 



-x 2 /2 



l-2e 1 



' 2'' 



; x ) + 2xv 



i-Pi 



l + 2e 1 
~~ 4~' 2 



r(^) 

r(3=2e) 



2e 3 



3 + 2e 3 



-; —x 
2' 



(5.1) 



where iiq is the confluent hyper geometric (Kummer) function. Note that, for e < 1/2, this 
solution will be nodeless for \u\ < 1 while it will have one node for \v\ > 1. 

Let us perform now a non-singular fc-th order SUSY transformation which creates precisely k 
new levels, additional to the standard ones E n = n + 1/2, n = 0, 1, 2, . . . of Hq, in the way 



Sp(H k 



1 3 

^ki • • • i ^i i Q ' ^ ' 



(5.2) 



where efc < • • • < ei < 1/2. In order that the Wronskian W(u\, . . . ,Uk) would be nodeless, the 
parameters i>i have to be taken as \v<\ < 1 for i odd and > 1 for i even, z = 1, . . . , k. The 
corresponding potential turns out to be 

V k (x) = ^-{ln[W( Ul ,...,u k )]}". 



It is important to note that there is a pair of natural ladder operators for H k : 

L k = B k a±B k' 
which are differential operators of (2k + l)-th order such that 

[H k ,L£\=±Li, 

and a~, a + are the standard annihilation and creation operators of the harmonic oscillator. 

From the intertwining relations of equation (|2.5p and the factorizations in equation f|2.4j) it 
is straightforward to show the following relation: 



Q2k+i(H k ) = L+L- = (h h - I) (H k - (H k - e t - 1 



i=i 



This means that {H k ,L k ,Lf} generate a (2fc)-th order PHA such that 
[L- k ,L + k ] = P 2k {H k ). 

Since the roots of Q2k+i(%) are {e k , . . . ,e±,l/2,e k + 1, . . . , e\ + 1}, it turns out that Sp(H k ) 
contains k one-step ladders, starting and ending at €j, j = 1, . . . , k, plus an infinite ladder which 
departs from 1/2 (compare with equation ()5.2[) ). 

Note that, for k = 1 the differential ladder operators L 1 = B^a^B^ are of third order, 
generating thus a second-order PHA. On the other hand, for k = 2 it turns out that = 
B+a ± B^ are of fifth order, and they generate a fourth-order PHA, etc. 
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6 Solutions of the Painleve IV equation through SUSY QM 

It was pointed out at the end of Section 0] the way in which solutions of Piv can be generated. 
Let us now employ this proposal for generating several families of such solutions. 

6.1 First-order SUSY QM 

We saw that, for k = 1, the ladder operators L,f are of third order. This means that the first- 
order SUSY transformation applied to the oscillator could provide solutions to the Piv equation. 
To find them, first we need to identify the extremal states, which are annihilated by and 
at the same time are eigenstates of H±. From the corresponding spectrum, one realizes that 
the transformed ground state of Hq and the eigenstate of H\ associated to t\ are two physical 
extremal states associated to our system. Since the third root of Q${x) is ei + 1 Sp(iii), then 
the corresponding extremal state will be nonphysical, which can be simply constructed from the 
nonphysical seed solution used to implement the transformation in the way Afa + u\. Due to 
this, the three extremal states for our system and their corresponding factorization energies (see 
equation (|4.ip ) become 

1 2 

ip £l oc — , ipe 2 °c Afe~ x /2 , ip£ 3 (x Afa + u 1: 

and 

Si = ex, £ 2 = -, £ 3 = ei + 1. 

The first-order SUSY partner potential V\(x) of the oscillator and the corresponding non-singular 
solution of Piv are 

9 

Vi(x) = X - - {hiM*)]}", g 1 (x,e 1 ) = -x- {ln[V> £l (*)]}' = -x + {ln[m(x)]}', (6.1) 

where we label the Piv solution with an index characterizing the order of the transformation 
employed and we indicate explicitly the dependence on the factorization energy. Notice that two 
additional solutions of the Piv can be obtained by cyclic permutations of the indices (1,2,3). 
However, they will have singularities at some points and thus we drop them in this approach. 
An illustration of the first-order SUSY partner potentials V\{x) of the oscillator as well as 
the corresponding solutions gi(x, e\) of Piv for e± = {l/4(blue), — 3/4(purple), — 7/4(red)} and 
v\ = 1/2 are shown in Fig. [2j 




-4 -2 2 4 -4 -2 2 4 



x x 

Figure 2. First-order SUSY partner potentials V\(x) (left) of the oscillator and the Piv solutions g\ (x, ei) 
(right) for ei = {l/4(blue), -3/4(purple), -7/4(red)} and v x = 1/2. 
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On the other hand, for k > 1 it is not clear that we can generate solutions of P\y. In the 
next subsection it will be shown that this is possible by imposing certain conditions on the seed 
solutions used to implement the SUSY transformation. 

6.2 fc-th order SUSY QM 

In Section [5] we saw that the k-th. order SUSY partners of the harmonic oscillator are ruled by 
(2&)-th order PHA. It is important to know if there is a way for this algebraic structure to be 
reduced to a second-order PHA and, if so, which are the requirements. The answer is contained 
in the following theorem. 

Theorem. Suppose that the k-th order SUSY partner H k of the harmonic oscillator Hamiltonian 
Hq is generated by k Schrddinger seed solutions Uj, j = l,...,k which are connected by the 
standard annihilation operator in the way: 

Uj = (a~Y~ l ui, ej = ei - (j - 1), j = l,...,k, (6.2) 

where u±(x) is a nodeless Schrddinger seed solution given by equation (|5.ip for ey < 1/2 and 

Wi\ < i. 

Therefore, the natural ladder operator = B^a + B^ of H k , which is of (2k + l)-th order, 
is factorized in the form 

L\ = P k ^(H k )l+, (6.3) 

where P k -\(H k ) = (H k — e±) ■ ■ ■ (H^ — e^-i) is a polynomial of (k — l)-th order in H^, it is 
a third-order differential ladder operator such that [Hk,lt] = it and 

lp- k = (H k - e k ) (n k - (H k -ex- 1). (6.4) 

Proof (by induction). For k = 1 the result is obvious since 
L+ = Po(H!)lf, P (#i) = 1. 

Let us suppose now that the theorem is valid for a given k; then, we are going to show that 
it is as well valid for k + 1. From the intertwining technique it is clear that we can go from H k 
to H k+ \ and vice versa through a first-order SUSY transformation 

H k +iA~l +l = A+ +1 H k , H k A^ +1 = A~ +1 H k . 

Moreover, it is straightforward to show that 

By using now the induction hypothesis of equation (|6.3[) for the index k it turns out that 

L k+i = A k+i P k-i( H k)l k A k+i = p k-i( H k+i)A^ +1 l^A^ +1 = P k -i(H k+1 )£+ +1 , (6.5) 

where £ k+ i = A'£, 1 l~£A'j~ +l is a fifth-order differential ladder operator for H k+ \. It is straight- 
forward to show 

^t+i^fc+i = ( H k+i ~ £k) 2 (H k+ i - e k+1 ) ^H k+1 - (H k+1 - e x - 1). 
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Note that the term (H k+ i — efc+i) [H k+ \ — |) (H k+ i — ei — 1) in this equation is precisely the 
result that would be obtained from the product ^ or ^ e third-order ladder operators 

of H k+ \. Thus, it is concluded that 

where q{H k+ \) is a polynomial of H k+ \. By remembering that i k+1 , l k+1 and -fffc+i are 
differential operators of 5-th, 3-th, and 2-th order respectively, one can conclude that q(H k +i) 
is lineal in H k -\-\ and therefore 

By substituting this result in equation (16.5|) we finally obtain: 

Lt +1 = P k -i(H k+1 )(H k+1 - e k )l+ +1 = P k (H k+1 )l+ +v U 

We have determined the restrictions on the Schrodinger seed solutions Uj to reduce the order 
of the natural algebraic structure of the Hamiltonian H k from 2k to 2. Now suppose we stick 
to these constraints for generating H k . Since the reduced ladder operator it is of third order, 
it turns out that we can once again obtain solutions of the Painleve IV equation. To get them, 
we need to identify the extremal states of our system. Since the roots of the polynomial of 
equation (j6.4|) are e k = e\ — (k — 1), 1/2, e\ + 1, the spectrum of H k consists of two physical 
ladders: a finite one starting from e k and ending at ei; an infinite one departing from 1/2. Thus, 
the two physical extremal states correspond to the eigenstate of H k associated to e k and to the 
mapped eigenstate of Hq with eigenvalue 1/2. The third extremal state (which corresponds to 
€\ + 1 Sp(-fffc)) is nonphysical, proportional to B^a + u\. Thus, the three extremal states are 

K W wu, ' ' ' ' ' ^ « B k e ~ x2/2 > ^3 « B+a + Ul , 
W\u\, ... ,u k ) 

with 

£ 1 = e k = e 1 -(k- 1), £ 2 = -, £3 = e\ + 1. 

The fc-th order SUSY partner of the oscillator potential and the corresponding non-singular 
solution of the Prv equation become 



2 

V k {x) = ^--{ln[W( Ul ,...,u k )}}", 
9k(x,ei) = -x - {ln[ip £l (x)]}' = -x - < In 



W{ui,. . . ,u k -{) 
W(ui,...,u k ) 



k>2. (6.6) 



Let us remind that the k Schrodinger seed solutions in the previous expressions are not longer 
arbitrary; they have to obey the restrictions imposed by our theorem (see equation (|6.2p ). 

We have illustrated the k-th order SUSY partner potentials V k {x) of the oscillator and the 
corresponding Pjy solutions g k (x, e±) in Fig. [3] (k = 2) and Fig. [H (k = 3) for ei = {l/4(blue), 
-3/4(purple),-7/4(red)} and v x = 1/2. 



7 Solution hierarchies. Explicit formulas 

The solutions g k (x,ei) of the Painleve IV equation can be classified according to the explicit 
functions of which they depend on [27]. Our general formulas, given by equations (|6.ip . (|6.6p . 
in general are expressed in terms of the confluent hypergeometric function \F±, although for 



10 



D. Bermudez and D.J. Fernandez C. 




Figure 3. Second-order SUSY partner potentials V^cc) (left) of the oscillator and the corresponding Prv 
solutions g2(x,ei) (right) for a = {l/4(blue), — 3/4(purple), — 7/4(red)} and v\ = 1/2. 




-4 -2 2 4 -4 -2 2 4 



x x 

Figure 4. Third-order SUSY partner potentials Vz(x) (left) of the oscillator and the corresponding Prv 
solutions gz(x,e\) (right) for ei = {l/4(blue), — 3/4(purple), — 7/4(red)} and v\ = 1/2. 

specific values of the parameter ei they can be simplified to the error function erf(x). Moreover, 
for particular parameters ei and they simplify further to rational solutions. 

Let us remark that, in this paper, we are interested in non-singular SUSY partner poten- 
tials and the corresponding non-singular solutions of -Prv- To accomplish this, we restrict the 
parameters to ei < 1/2 and \v\\ < 1. 



7.1 Confluent hypergeometric function hierarchy 



In general, the solutions of Pry given in equations (|6.1[) . ()6.6[) are expressed in terms of two 
confluent hypergeometric functions. Next we write down the explicit formula for gi(x,€\) in 
terms of the parameters e\, v\ (with e\ < 1/2 and \i>i\ < 1 to avoid singularities): 

2utT [(3 - 6x 2 ) |; x 2 ) + x 2 (3 - 2 Cl ) 1 F 1 (Z^a, §; x 2 )] 



3r (1=*L) lFl (l^a, i; x 2 ) + 6^xr (^i) 1 P 1 (3^1, |;x 2 ) 
3xF (1^) [-2 1 P 1 (±^a I; x 2 ) + (1 - 2 £1 ) 1 P 1 («=*i, |; x 2 )] (7 i) 



3r (!=*!) 1 f i (1=*l, i; x 2 ) + e^xr (2^a) 1 p 1 (2^a, |. x2) • 

Notice that the solutions g\{x,e\) plotted in Fig. [2] for specific values of ei, i^i correspond 
to particular cases of this hierarchy. The explicit analytic formulas for higher-order solutions 
<7fc(x, ei) can be obtained through formula (|6.6p . and they have a similar form as in equation (|7.ip . 
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Figure 5. First-order SUSY partner potentials V\(x) (left) of the oscillator and the Pry solutions 
gi(x, ex) (right) for e\ = {—1/2 (blue), —3/2 (purple), —5/2 (red)} and v\ = 0.999, which belong to the 
error function hierarchy of solutions. 



7.2 Error function hierarchy 

It would be interesting to analyze the possibility of reducing the explicit form of the Pry solution 
to the error function. To do that, let us fix the factorization energy in such a way that any of 
the two hypergeometric series of equation (15.ip reduces to that function. This can be achieved 
for ex taking a value in the set 



(2m + 1) 



(7.2) 



If we define the auxiliary function (p ui (x) = ^pne x [1 + v\ erf(x)] to simplify the formulas, we 
can get simple expressions for gk(x,ei) with some specific parameters k and e\: 





-1/2) = 


9l(x, 


-3/2) = 


9l(x, 


-5/2) = 


92(x, 


-1/2) = 



2V! 

¥Vi 0) ' 

1 + xip Ul (x) ' 

4hvi + <p„ 1 (a)] 
2v l x + {l + 2x 2 )vv l 0)' 

4z/i[^i + 6ip Ul (x)] 

2uf) 



(7.3) 



(7.4) 



<Pm (x) (x) - 2v lX ip Ul (x) 



An illustration of the first-order SUSY partner potentials V\{x) of the oscillator and the 
corresponding Piv solutions gi{x,e\) of equations (j7.3|) . (|7.4p is given in Fig. [5j 



7.3 Rational hierarchy 

Our previous formalism allows us to generate solutions of Piv involving in general the confluent 
hypergeometric series, which has an infinite sum of terms. Let us look for the restrictions needed 
to reduce the explicit form of equation (|6.6p to non-singular rational solutions. To achieve this, 
once again the factorization energy e\ has to take a value in the set given by equation (17. 2p . 
but depending on the ei taken, just one of the two hypergeometric functions is reduced to 
a polynomial. Thus, we need to choose additionally the parameter v\ = or v\ — > oo to keep 
the appropriate hypergeometric function. However, for the values —(4m — l)/2, m = 1,2,... 
and v\ — > oo, it turns out that u\ will have always a node at x = 0, which will produce one 
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Figure 6. Third-order SUSY partner potentials Va(x) (left) of the oscillator and the Piy solutions 
g 3 (^,ei) (right) for e x = {-5/2(blue), -9/2(purple), -13/2(red)} and v x = 0. 



singularity for the corresponding Pry solution. In conclusion, the rational non-singular solutions 
gk(x, ei) of the Piv arise by making in equation ()5.ip v\ = and 



ei € 



(4m + 1) 



Taking as the point of departure Schrddinger solutions with these v\ and ei and using our 
previous expressions (|6.6p for a given order of the transformation we get the following explicit 
expressions for gi-(x, ei), some of which are illustrated in Fig. [6j 



Si(x,-5/2) = 
<7i Or, -9/2) = 
5i (x, -13/2) 
S2(x,-5/2) = 
g 2 {x,-9/2) = 
g 2 (x,-13/2) 

S3 (x, -5/2) = 
S3 (x, -9/2) = 



4x 



1 + 2x 2 ' 

8(3x + 2x 3 ) 
3 + 12x 2 + 4x 4 ' 

12(15a + 20x 3 + 4x 5 ) 
" 15 + 90x 2 + 60x 4 + 8x 6 ' 

4x 16x 3 
~ 1 + 2x 2 + 3 + 4x 4 ' 

8(3x + 2x 3 ) 32(15x 3 + 12x 5 + 4x 7 ) 

~3 + 12x 2 + 4x 4 + 45 + 120x 4 + 64x 6 + 16x 8 ' 

12(15x + 20x 3 + 4x 5 ) 
" ~ 15 + 90x 2 + 60x 4 + 8x 6 

48(525x 3 + 840x 5 + 600x 7 + 160x 9 + 16X 11 ) 
1575 + 6300x 4 + 6720x 6 + 3600x 8 + 768x 10 + 64x 12 ' 
4x(27 - 72x 2 + 16x 8 ) 



+ 



27 + 54x 2 + 96x 6 - 48x 8 + 32x 10 ' 

32(15x 3 + 12x 5 + 4x 7 ) 
~ 45 + 120x 4 + 64x 6 + 16x 8 

24(225x - 150x 3 + 120x 5 + 240x 7 + 80x 9 + 32X 11 ) 
675 + 2700x 2 - 900x 4 + 480x 6 + 720x 8 + 192x 10 + 64x 12 ' 



+ 



8 Conclusions 



In the first part of this paper we have reviewed the main results concerning the most gen- 
eral Schrddinger Hamiltonians characterized by second-order PHA, i.e., possessing third-order 
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differential ladder operators. In particular, it was seen that the corresponding potentials can be 
obtained from the solutions to P\y. 

On the other hand, starting from the k-th order SUSY partners of the harmonic oscillator 
potential, a prescription for generating solutions of Ppv has been introduced. We have shown 
that the Hamiltonians associated to these solutions have two independent physical ladders: an 
infinite one starting from 1/2 and a finite one placed completely below 1/2. We also have 
identified three solution hierarchies of the P\y equation, namely, confluent hyper geometric, 
error function, and rational hierarchies, as well as some explicit expressions for each of them. 

Inside the idea of spectral manipulation, it would be interesting to investigate the possibility 
of constructing potentials with more freedom for the position of the finite physical ladder, 
e.g., some or all levels of this ladder could be placed above 1/2. This is a subject of further 
investigation which we hope to address in the near future. 
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